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Kibble and Zurek have provided a unifying picture for the 
onset of phase transitions in relativistic QFT and condensed 
matter systems respectively, strongly supported by agreement 
with condensed matter experiments in 3 He. The failure of 
a recent experiment on 4 He to agree with the predictions 
has prompted a reappraisal of this picture. We provide an 
alternative picture that explains the experimental evidence. 



I. THE ONSET OF PHASE TRANSITIONS 

The two environments in which thermal field theory is 
relevant are the early universe and heavy-ion collisions, 
in both of which there are changes of phase. Not surpris- 
ingly, one of the themes of this series of Workshops on 
Thermal Field Theory (TFT) has been phase transitions. 

Consider a situation in which the symmetry group of 
the theory is broken, on cooling through the critical tem- 
perature T c , by its degenerate groundstate manifold. Us- 
ing the tools of equilibrium TFT we can determine the 
nature of the transition. At late times after the transition 
the fields are ordered on large scales, in that the field will 
adopt a single value from this degenerate set over a large 
spatial region. This might seem to be all that we need 
to know, but how this is achieved can have interesting 
consequences. 

Such problems require that we go beyond equilibrium 
TFT. In practice, we often know remarkably little about 
the dynamics of thermal systems. For simplicity, I shall 
assume scalar field order parameters, with continuous 
transitions. In particular, I want to discuss some aspects 
of the onset of such phase transitions, the very early times 
after the implementation of a transition when the scalar 
fields are only just beginning to become ordered. 

A simple question to ask is the following: In princi- 
ple, the field correlation length diverges at a continu- 
ous transition. In practice, it does not. What happens? 
This is relevant for transitions that leave topological de- 
fects like walls, monopoles, vortices, or textures in their 
wake since we might expect 'defects' to be just that, enti- 
ties whose separation is characterised by the correlation 



length. If this were simply so, a measurement of defect 
densities would be a measurement of correlation lengths. 
Estimates of this early field ordering have been made by 
Kibble [fj]j^], using simple causal arguments, because of 
the implication of defects for astrophysics. Vortices, in 
particular (cosmic strings), can be important for struc- 
ture formation in the universe. [|] 

There are great difficulties in converting predictions for 
the early universe into experimental observations. Zurek 
suggested |J that similar arguments were applicable to 
condensed matter systems for which direct experiments 
on defect densities could be performed. This has lead to a 
burst of activity from theorists working on the boundary 
between QFT and Condensed Matter theory and from 
condensed matter experimentalists. To date almost all 
experiments have involved superfluids, with scalar order 
parameters of the type considered above and all but one 
experiment is in agreement with these simple causal pre- 
dictions. 

In this talk I shall 

• review the Kibble/Zurek causality predictions for 
initial correlation lengths and defect densities and 
summarise the results of the condensed matter ex- 
periments. 

• present an alternative picture for the onset of defect 
production, applicable to both relativistic quantum 
fields and condensed matter systems. 

• show how this alternative picture gives essentially 
the same results as the Zurek picture for those con- 
densed matter systems for which there is experi- 
mental agreement. 

• provide an explanation for why some condensed 
matter experiments will be in disagreement with 
Zurek's predictions, including the experiment in 
question. 

• try to draw some broader conclusions about the rel- 
evance of Kibble's predictions for phase transitions 
in QFT. 



*Talk presented at 5th International Workshop on Thermal 
Field Theories and Their Applications, Regensberg, August 
1998 



x In a different context, the topological textures of the NLSM 
0(4) theory can be identified with Skyrmionic hadrons pro- 
duced in the hadronisation of quark-gluon plasma. 
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II. WHEN SYMMETRY BREAKS, HOW BIG 
ARE THE SMALLEST IDENTIFIABLE PIECES? 

Insofar that defect density can be correlated simply to 
correlation length, the maximum density (an experimen- 
tal observable in condensed matter systems, although not 
for the early universe) will be associated with the smallest 
identifiable correlation length in the broken phase once 
the transition has been effected. This provides the initial 
condition for the evolution of field ordering. In order to 
see how to identify these 'smallest pieces]^] it is sufficient 
to consider the simplest theory, that of a single relativis- 
tic quantum scalar field in three spatial dimensions with a 
double well potential, undergoing a temperature quench. 
In the first instance we assume that the qualitative dy- 
namics of the transition are conditioned by the field's 
equilibrium free energy, of the form 

F(T) = J d 3 x Q(V</>) 2 + \m 2 {T)d? 2 + \\A . (1) 

Prior to the transition, at temperature T > T c , the crit- 
ical temperature, m(T) > plays the role of an effective 
'plasma' mass due to the interactions of <p with the heat- 
bath, which includes its own particles. After the transi- 
tion, when T is effectively zero, m 2 (0) = — M 2 < en- 
forces the Z2 symmetry-breaking, with field expectation 
values (4>) — ±77, rj 2 = M 2 /X. The change in temperature 
that leads to the change in the sign of m 2 is most simply 
understood as a consequence of the system expanding. 
Thus, in the early universe, a weakly interacting rela- 
tivistic plasma at temperature T > M has an entropy 
density s oc T 3 . As long as thermal equilibrium can be 
maintained, constant entropy S per comoving volume, 
S oc sa(t) 3 , gives T oc a(t) -1 and falling, for increasing 
scale factor a(t). We shall be too simple in our discussion 
to warrant the inclusion of a metric in Eq.|l|. 

Some caution is needed in interpreting V((j>) = 
m 2 (T)4> 2 /2 + X(j} 4 /4: as the conventional Effective Po- 
tential since this is a zero-momentum, or infinite- volume 
field average, construct. At the very least, a more realis- 
tic potential for dynamics would be coarse-grained to the 
relevant scale. In practice, as long as we never look on 
scales smaller that the field correlation length this turns 
out not to matter much. The only relevant point is that 
m 2 (T) vanishes at T — T c , true however V((f>) is defined. 

The purpose of V(<j>, T) is, initially, to set the scale 
of the fluctuations of cj> about = 0, the subsequent 
false vacuum. The equilibrium correlation length of the 
field fluctuations is £ eg (T) = |m _1 (T)|. It is sufficient 
to adopt a mean-field approximation in which m 2 (T) oc 
(T—T c ). For fields averaged on the scale of £ e q(T) we find 
(4> 2 ) = 0(m(T)T), much less than rj 2 for our relativistic 



2 The title of this section is essentially that posed in recent 
papers by Zurek Q. 



theory in which T c = 0(rj). After the transition the field, 
at each point in space, begins to collapse to one of the 
false vacua cf> = ±77, rj 2 = Af 2 /A, as signalled by V(4>), 
initially in a random way. However, the field soon begins 
to become organised into 'domains', in each of which the 
field has a constant sign. The defects of this simple model 
are domain walls, across which the field flips from one 
minimum to the other. Defect density then looks to be 
related to domain size. 

How this collapse takes place determines the size of 
the first identifiable domains. It was suggested by Kibble 
and Zurek that this size is essentially the field correlation 
length £ eg at some appropriate temperature close to the 
transition. I shall argue later that this is too simple but, 
nonetheless it is a plausible starting point. Near the con- 
tinuous transition £ eq becomes large, approaching infinity 
at T c . In practice physical correlations remain finite. The 
question then becomes: how large does the field correla- 
tion get in practice? Two very different mechanisms have 
been proposed for estimating this size. 

A. Thermal activation 

In the early work on the cosmic string scenario it 
was assumed 0] that initial domain size was fixed in 
the Ginzburg regime. By this we mean the following. 
Suppose the temperature T(t) varies sufficiently slowly 
with time t that it makes sense to replace V(cf>,T) by 
V{4>, T(t)). Well away from the transition this is justified, 
but close to the transition it is not. A band in tempera- 
ture about T c in which it is not sensible is T G < T < T G , 
where T G are the Ginzburg temperatures. As we cool 
down in the vicinity of T c , V(<fi,T) ceases to be a re- 
liable guide to the field fluctuations at T G , where the 
relative Ginzburg temperature is of order of magnitude 
|1 — /T c \ = 0(A). Once we are below T c , and the cen- 
tral hump in V(4>,T(t)) is forming, T G signals the tem- 
perature above which there is a significant probability for 
thermal fluctuations over the central hump on the scale 
of the correlation length. Most simply, it is determined 
by the condition 

^V{T G )e eq {T G )^T G (2) 

where AV(T) is the difference between the central max- 
imum and the minima of V(cj>,T). Again we find 
\1-T G /T C \=0(\). _ 

Whereas, above (T G ) there will be a population of 'do- 
mains', fluctuating in and out of existence, at tempera- 
tures below T G fluctuations from one minimum to the 
other become increasingly unlikely. When this happens 
the correlation length is 

= 000, (3) 

where £0 = Af _1 is the natural unit of length, the Comp- 
ton wavelength of the <f> particles. 
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It is tempting to identify £ eg (Tg) with the scale at 
which stable domains begin to form. Subsequent work 
shows this to be incorrect. Thermal fluctuations are rel- 
evant to the formation of small domains, and to wiggles 
in the boundaries of larger domains, but not in the for- 
mation of larger domains themselves. That is an issue 
that requires more than equilibrium physics. The most 
simple dynamical arguments can be understood in terms 
of causality. 



B. Causality 

If the Ginzburg criteria attempt to set scales once the 
critical temperature has been passed, causal arguments 
attempt to set scales before it is reached. We have seen 
that £ eq (T(i)) diverges at Tit) — T c , which we suppose 
happens at t = 0. This cannot be the case for the true 
correlation length £(t), which can only grow so far in a 
finite time. Initially, for t < 0, when we are far from 
the transition, we again assume effective equilibrium, 
and the field correlation length £(t) tracks £ eq (T(t)) ap- 
proximately. However, as we get closer to the transition 
£ e g(T(t)) begins to increase arbitrarily fast. As a crude 
upper bound, the true correlation length fails to keep up 
with £ eg (T(i)) by the time — tc at which £ e(J is growing 
at the speed of light, d£ eq (T(—tc))/dt = 1. It was sug- 
gested by Kibble Q that, once we have reached this time 
£(t) freezes in, remaining approximately constant until 
the time t ~ +tc after the transition when is once again 
becomes comparable to the now decreasing value of £ eq . 
The correlation length £, eq (tc) = £e<?(— tc) is argued to 
provide the scale for the minimum domain size after the 
transition. 

Specifically, if we assume a time-dependence m 2 (t) = 
—M 2 t/tQ in the vicinity of t — 0, when the transition 
begins to be effected, then the causality condition gives 



t c = £q /3 (2M;t 2 / 3 . As a result, 



M^ eq {t c ) = (Mr ) 1/3 , 



which, with condensed matter in mind, we write as 

- > V3 
\ T 



(4) 



(5) 



where tq — £o = M^ 1 are the natural time and distance 
scales. In contrast to Eq|], Eq.||depends explicitly on the 
quench rate, as we would expect. For tq 3> tq the field 
is correlated on a scale of many Compton wavelengths. 



C. QFT or Condensed Matter 

This approach of Kibble was one of the motivations 
for a similar analysis by Zurek || of transitions with 
scalar order parameters in condensed matter. Qualita- 
tively, neither the Ginzburg thermal fluctuations, with 



fluctuation length Eq.||, nor the simple causal argument 
above depend critically on the fact that the free energy 
Eq.|l| is originally assumed to be derived from a relativis- 
tic quantum field theory. After rescaling, F could equally 
well be the Ginzburg-Landau free energy 



F(T) 



d 6 x 



2m 



(V0) 2 + a(T)</> 2 



(6) 



for a non-relativistic condensed matter field, in which 
a(T) ex m 2 {T) vanishes at the critical temperature T c . 
The only difference is that, in the causal argument, the 
speed of light should be replaced by the speed of sound, 
with different critical index. 

Explicitly, let us assume the mean-field result a(T) = 
aQe(T c ), where e = (T/T c — 1), remains valid as T/T c 
varies with time t. In particular, we first take a(t) = 
a(T(t)) = —aot/rQ in the vicinity of T c . Then the fun- 
damental length and time scales £o and to are given from 
Eq.^as £q = h 2 /2mao and tq = h/ctQ. It follows that the 
equilibrium correlation length £ e9 (i) and the relaxation 
time r(t) diverge when t vanishes as 



£e 9 (*)=£c 



t 



-1/2 



r(t) = to 



f 



(7) 



The speed of sound is c(t) = t; eq (t)/T(t), slowing down as 
we approach the transition as It] 1 / 2 . The causal coun- 
terpart to d£ eq (t)/dt — 1 for the relativistic field is 
d£eq(t)/dt — c(t). This is satisfied at t = —tc, where 
tc = y/TQTQ, with corresponding correlation length 

/ \ 1/4 

£e 9 (*c) = £e 8 (-M=£of^J • (8) 

(cf. Eq.^|). A variant of this argument that gives es- 
sentially the same results is obtained by comparing the 
quench rate directly to the relaxation rate of the field 
fluctuations. We stress that, yet again, the assumption 
is that the length scale that determines the initial corre- 
lation length of the field freezes in before the transition 
begins. Whatever, the field is already correlated on a 
scale of many Compton wavelengths when it begins to 
unfreeze. 



III. EXPERIMENTS 

The end result of the simple causality arguments is 
that, both for QFT and condensed matter, when the field 
begins to order itself its correlation length has the form 



£e 9 (ic)=M — 



(9) 



for appropriate 7. ^ 



3 In fact, the powers of Eq.^ and EqJ^ are mean-field results, 
changed on implementing the renormalisation group. 
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For a relativistic complex scalar field <p with a global 
U(l) or 0(2) symmetry the causality argument is the 
same, but the physical situation is very different. Con- 
sider the theory controlled by a free energy 

F(T) = J d 3 x (W| 2 + m 2 (T)|0| 2 + A|0| 4 ) (10) 

with m 2 (T) switching from positive to negative. The 
minima of the final potential of Eq.|l^ now constitute the 
circle <fi = rje 10 " . When the transition begins <fi begins to 
fall into the valley of the potential, choosing a random 
phase. This randomly chosen phase will vary from point 
to point. Solutions to 8F/8<fi — include vortices, topo- 
logical defects, tubes of 'false' vacuum <f> « 0, around 
which the field phase changes by ±2ir. In an early uni- 
verse context these are 'cosmic strings', but we shall not 
consider their properties here. The jump that Kibble 
made was to assume that the correlation length EqJ|, 
equally applicable to a complex field, also sets the scale 
for the typical minimum intervortex distance. 
That is, the initial vortex density ridef is 

-'=°(w)=°(3ter> <n) 

for 7 — 1/3. Equivalently, the length of vortices in a box 
volume v is 0(n<i e fv). We stress that this assumption is 
independent of the argument that lead to Eq.[5]. 

Since £o also measures cold vortex thickness, tq ^> tq 
corresponds to a measurably large number of widely sep- 
arated vortices. If it could be argued that this initial 
network behaves classically then, thereafter, the density 
will reduce due to the collapse of small loops, intersec- 
tions chopping off loops which in turn collapse, and vor- 
tex straightening so as to reduce the gradient energy of 
the field. 

Even if cosmic strings were produced in so simple a way 
in the very early universe it is not possible to compare the 
density Eq.|ll| with experiment. It was Zurek who first 
suggested that this causal argument for defect density be 
tested in condensed matter systems, particularly in liquid 
helium. 



A. Superfluid helium 

Vortex lines in both superfluid He and 3 He are good 
analogues of global cosmic strings. A crude but effective 
model is to treat the system as composed of two fluids, 
the normal fluid and the superfluid, which has zero vis- 
cosity. In A He the bose superfluid is characterised by a 
complex field 0, whose squared modulus \4>\ 2 is the su- 
perfluid density. The superfluid fraction is unity at abso- 
lute zero, falling to zero as the temperature rises to the 
lambda point at 2.17K. The Landau- Ginzburg theory for 
A He has, as its free energy, 



F(T) = Jd 3 x (^|l|V0| 2 +a(T)|0| 2 + i/3|^, (12) 

the scaled counterpart of Eq.[l(]. The static classical field 
equation SF/S(j) — has vortex solutions as before, with 
width £o- 

The situation is more complicated, but more inter- 
esting, for 3 He, which becomes superfluid at the much 
lower temperature of 2 mK. The reason is that the 3 He 
is a fermion. Thus the mechanism for superfluidity is 
very different from that of i He. Somewhat as in a BCS 
superconductor, these fermions form the counterpart to 
Cooper pairs. However, whereas the (electron) Cooper 
pairs in a superconductor form a 1 S state, the 3 He pairs 
form a 3 P state. The order parameter A a i is a com- 
plex 3x3 matrix A a i. There are two distinct superfluid 
phases, depending on how the 50(3) x SO(3) x U(l) 
symmetry is broken. If the normal fluid is cooled at low 
pressures, it makes a transition to the 3 He — B phase, in 
which A a i takes the form A a i = i? Q i(oj)e 1 *, where R is a 
real rotation matrix, corresponding to a rotation through 
an arbitrary u> J5| 

The Landau-Ginzburg free energy is, necessarily, more 
complicated, but the effective potential V(A ai , T) has the 
diagonal form V(A,T) = a(T)\A ai \ 2 + 0(A 4 ) for small 
fluctuations, and this is all that we need for the pro- 
duction of vortices at very early times. Thus the Zurek 
analysis leads to the prediction Eqjll], as before, for ap- 
propriate 7. However, for 3 He the mean-field approxi- 
mation is good and the mean- field critical index 7 = 1/4 
is not renormalised, whereas for A He a better value is 
7 = 1/3. 

B. Experiments in 3 He. 

Although 3 He is more complicated to work with, the 
experiments to check Eq.Jll] are cleaner for both exper- 
imental and theoretical reasons. First, because the nu- 
cleus has spin 1/2, even individual vortices can be de- 
tected by magnetic resonance. Second, because vortex 
width is many atomic spacings the Landau-Ginzburg the- 
ory is good. 

So far, experiments have been of two types. In the 
Helsinki experiment Q superfluid 3 He in a rotating cryo- 
stat is bombarded by slow neutrons. Each neutron en- 
tering the chamber releases 760 keV, via the reaction 
n + 3 He — > p + 3 He + 760keV. The energy goes into 
the kinetic energy of the proton and triton, and is dissi- 
pated by ionisation, heating a region of the sample above 
its transition temperature. The heated region then cools 
back through the transition temperature, creating vor- 
tices. Vortices above a critical size (dependent on the 
angular velocity of the cryostat) grow and migrate to the 
centre of the apparatus, where they are counted by an 
NMR absorption measurement. Suffice to say that the 
quench is very fast, with tq/tq — O(10 3 ). Agreement 
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with Eq.|Tl| and Eq.^ is very good, at the level of an or- 
der of magnitude. 

The second type of experiment has been performed 
at Grenoble and Lancaster Jjj. Rather than count indi- 
vidual vortices, the experiment detects the total energy 
going into vortex formation. As before, 3 He is irradiated 
by neutrons. After each absorption the energy released 
in the form of quasiparticles is measured, and found to 
be less than the total 760 keV. This missing energy is 
assumed to have been expended on vortex production. 
Again, agreement with Zurek's prediction Eq.Jll] and Eq.|8| 
is good. 



C. Experiments in He. 

The experiments in 4 iJe, conducted at Lancaster, fol- 
low Zurek's original suggestion. The idea is to expand a 
sample of normal fluid helium, in a container with bel- 
lows, so that it becomes superfluid at essentially constant 
temperature. That is, we change 1 — T/T c from negative 
to positive by reducing the pressure, thereby increasing 
T c . As the system goes into the superfluid phase a tangle 
of vortices is formed, because of the random distribution 
of field phases. The vortices are detected by measur- 
ing the attenuation of second sound within the bellows. 
Second sound scatters off vortices, and its attenuation 
gives a good measure of vortex density. A mechanical 
quench is slow, with tq some tens of milliseconds, and 
tq/tq — O(10 10 ). Two experiments have been performed 
P,P| . In the first fair agreement was found with the pre- 
diction Eq.|ll|, although it was not possible to vary tq. 
However, there were potential problems with hydrody- 
namic effects at the bellows, and at the capillary with 
which the bellows were filled. A second experiment, de- 
signed to minimise these and other problems has failed 
to see any vortices whatsoever. 

There is certainly no agreement, in this or any other 
experiment on 3 He, with the thermal fluctuation density 
that would be based on Eq.|[ 



IV. THE KIBBLE-ZUREK PICTURE FOR THE 
FREEZING IN OF £ IS CORRECT. 

We are in the strange position that, if the second Lan- 
caster experiment Q is correct, either both the Zurek 
predictions Eq.|ll] and EqJ| works well, or one or both 
fail. We shall now argue that the freezing in of £ as in 
Eq.|| is correct, essentially on dimensional grounds. To 
show this we need a concrete model for the dynamics. 



A. Condensed matter: the TDLG equation 

We assume that the dynamics of the transition can 
be derived from the explicitly time-dependent Landau- 



Ginzburg free energy 



a(W| 2 + i /?H 4 



(13) 



In (H) = (0i+i0a)/V2 is the complex order-parameter 
field, whose magnitude determines the superfluid density. 
In equilibrium at temperature T, in a mean field ap- 
proximation, the chemical potential a(T) takes the form 
a(T) — Qo e (7c), where e = (T/T c — 1). In a quench in 
which T c or T changes it is convenient to shift the origin 
in time, to write e as 



e(t) = e 



-9{t) 



(14) 



for — oo < t < tq(1 + eo), after which e(t) = —1. eo 
measures the original relative temperature and tq defines 
the quench rate. The quench begins at time t = and 
the transition from the normal to the superfluid phase 
begins at time t = €qtq. 

Motivated by Zurek's later numerical jl| simulations 
with the time-dependent Landau-Ginzburg (TDLG) 
equation for F of (|l^), we assume a linear response 

1 »> S " ' *, (15) 



r dt 



where ri a is Gaussian noise. More details are given else- 



where 1 10 . We can show self-consistently that, for the 
relevant time-interval — tc < At < tc the self-interaction 
term can be neglected (f3 = 0). This both preserves 
Gaussian field fluctuations and leads to £ « ^(tq/tq) 1 ^ 4 
arising in a natural way, as we shall see. 

At relative time At the diagonal equal-time correlation 



function G(r, At) is defined by 

{M r )M0))At = SdbG{r,At). 



(16) 



Anticipating that correlation lengths are approxi- 
mately frozen in during this period it is sufficient to per- 
form our calculations at At = 0. Eq.[l5| is solvable, in 
time and space units tq and £oj to give the Fourier trans- 
form of G(r, At) as 



G(k) 



dte- t2/TQ : 



T Q k 2 



(17) 



independent of eo, provided €qTq, c§tq ^ 1. For 4 He, 
where this is more problematical, this is valid, since al- 



io 



10 



is so 



though e ~ 10 — 10 is very small, tq 
large. 

On dimensional grounds the correlation length of the 
field is 0(^ q (t c )) = 0(eo(T Q /r ) 1 / 4 ). However, rather 
than have an asymptotic fall-off of the form e~ T l^" t ^ tc ^ 1 
for large r, G(r) oc exp(— 0((r/£ e(? (t c )) 4 / 3 )). Nonethe- 
less, numerically, it is remarkably well represented by 
e - r /?e< ! ( t c) ; w ith coefficient unity in the exponent, for r 
being a few multiples of £, e q{tc)i for reasons that are not 
clear to us. In that sense Zurek's prediction is robust, 
since explicit calculation j2l|] shows that £(t) does not 
vary strongly in the interval — tc < At < tc- 
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B. QFT: the free roll 

The dynamical equations for a hot quantum field 
(Heisenberg equations with Boltzmann boundary condi- 
tions) are very different from the empirical TDLG equa- 
tion above. Fortunately, as for the condensed matter 
case, the interval — tc < At < tc occurs in the linear 
regime, when the operator equations 



O 2 



SF 



or 2 



(18) 



for F of Eq|T(], are solvable in terms of the mode functions 
Xfc (*)) satisfying 



alt 2 



subject to Xk M 



m 2 (t) 



(19) 



at t < 0, for incident fre- 
quency uiin — Vk 2 + £qM 2 , for m 2 (t) = e(t)M 2 , where 
e(t) is parametrised as for the TDLG equation above. 
This corresponds to a temperature quench from an ini- 
tial state of thermal equilibrium at temperature T > T c , 
where (Tq/T c — 1) = eo- There is no reason to take eo 
small. The diagonal correlation function G(r, t) of Eq.|l6| 
is given as the equaltime propagator 



G(r,t) = 



_ik.x~ A- 



xt(t) Xk -(t)c(k), 



(20) 



where C(k) = coth(cJi„(fc) /2Tq) /2u>i n (k) encodes the ini- 
tial conditions. 

An exact solution can be given [ p| in terms of Airy, 
and related, functions. Dimensional analysis shows that, 
at displaced time At = 0, £ eq (tc) of EqJ^ again sets the 
scale of the equaltime correlation function. We have yet 
to work out whether the coefficient is appproximately 
unity, but Kibble's insight is correct. Further detail is 
unnecessary since we shall argue that it is the assumption 
that we can infer the defect density from £ e g(icO that is 
flawed. 



quench rate for 3 He experiments is several orders of mag- 
nitude faster than the mechanically-driven quenches of 
the 4 He experiments. Secondly, for 3 He the Ginzburg 
regime is extremely narrow, whereas for 4 He it is very 
broad. In fact, the A He experiments begin and end in 
the Ginzburg regime, where thermal fluctuations are im- 
portant. The causality arguments arc too simple to ac- 
comodate these facts. 

If these differences are to be visible in the formalism, 
it can only be through the way in which we relate vor- 
tex density to correlation length. That is, the flaw of 
the prediction lies in the assumption Eq.O, as becomes 
apparent once we appreciate that it hides the problem 
of how to count vortices. One way is to work through 
line zeroes. For convenience we restrict ourselves to 0(2) 
vortices whose cores are line zeroes of the complex field 
<$>. Q The converse is not true since zeroes occur on all 
scales. However, a starting-point for counting vortices 
in superfluids is to count line zeroes of an appropriately 
coarse-grained field, in which structure on a scale smaller 
than £o, the classical vortex size, is not present |0|. This 
is also the unstated basis of the numerous numerical sim- 
ulations p3| of cosmic string networks built from Gaus- 



sian fluctuations (but see |L4 



For the moment, we put 



in a cutoff I = 0(£o) by hand, as 



Gi(r,£t)= / fke lkx G(k,At)e 



(21) 



We stress that the long-distance correlation length 
£,eq(tc) depends essentially on the position of the nearest 
singularity of G(fc, At) in the complex k-plane, indepen- 
dent of I. 

This is not the case for the line-zero density n zero , 
depending, in our Gaussian approximation jl5 16 1, on the 
short- distance behaviour of G;(r, At) as 



n zero (At) 



-lG?(0,At) 



(22) 



2tt Gi(0,At) ' 

the ratio of fourth to second moments of G(k, At) e~ fe2/ 



V. DEFECT DENSITIES DO NOT DETERMINE 
CORRELATION LENGTHS DIRECTLY 

We have seen that there is no reason to disbelieve the 
causal arguments of Kibble for QFT and Zurek for con- 
densed matter as to the time at which the field fluctu- 
ations freeze in, and the associated correlation length. 
The excellent agreement with the 3 He experiments also 
shows that, for condensed matter, this length does, in- 
deed, characterise vortex separation at the time when the 
defects form. 

However, if we take the Lancaster experiment at face 
value, this cannot always be the case. If we think about 
why this might be, the differences between the 3 He and 
4 He experiments are twofold. Firstly, the nuclear-driven 



A. TDLG condensed matter 

In the linear regime everything is calculable. If we 
define the line-zero separation £ zero by 



n zero (At) = 



1 



27T^ ero (Ai) 2 



(23) 



it is apparent that ^ zer o(At) has little, if anything, to do 
with £ eq (Ai) directly. 



4 More complicated vortices require more complicated anal- 
ysis, but we expect essentially the same conclusions. 



G 



In fact, at At 



2 

zero 



(24) 



independent of eo- We have a situation in which the den- 
sity of line zeroes depends entirely on the scale at which 
we look. Such fractal behaviour cannot be understood 
as representing stable vortices, although such line zeroes 
might be termed fluctuation vortices. We would not wish 
to identify the length scale £ 2eT - with defect separation 
at this time, even if the latter could be defined. 

This is not surprising. Although the field correlation 
length £(At) may have frozen in by At — 0, the symme- 
try breaking has not been effected. In a classical sense at 
least, vortices can only be identified once the field mag- 
nitude has grown to its equilibrium value 



(H 2 > = ao/P, 



(25) 



and we should not begin to count them before then. Even 
though the long-range correlation length £(At) will not 
change substantially in that time, this will not be the case 
for £ zero (At), as long-range modes increase in amplitude 
as the field becomes ordered. 

Specifically, on continuing to use for times t > 
£q t q we see that, as the unfreezing occurs, long wave- 
length modes with k 2 < tjrQ — cq grow exponentially. 
Provided £§tq ^> 1 they soon begin to dominate the cor- 
relation functions. Let 



G n (At) 



dkk 2n G(k,At) 



(26) 



be the moments of G(k, At). In units of £o and r o we 
find PI that 



G n (At) 



n_ JAt/t} 2 



2™+!/2 



dt' 



g-(t'-At) 2 /* 2 
[V + l 2 /2] n + 1 / 2 ' 



(27) 



where we measure the dimensionless time At in units of 
t = y/TQ from t — e TQ and /„ = J Q dkk 2n e~ k . For 
small relative times the integrand gets a large contribu- 
tion from the ultraviolet cutoff dependent lower endpoint, 
and we recover (Q). As long as the endpoints make 
a significant contribution to the whole then the density 
of line zeroes derived from (|27|) will be strongly depen- 
dent on scale. Only when their contribution is small and 
dn zero /dl is small in comparison to n zero /l at I = £o 
can we identify the essentially non-fractal n zero with a 
meaningful vortex density. 

Ati and At2 , the times at which the exponential modes 
begin to dominate in the integrands of G\ and G2, can 
be determined by comparing the relative strengths of 
the contributions from the scale-independent saddlepoint 
and the endpoint in (p7|). If p\ — Ati/i and P2 = A^/f 
are the multiples of i at which it happens then, for the for- 



mer to dominate the latter requires e Pl /p 1 
' - 5/2 >3v^ 



3/2 > 



7V27T 



and e P2 /p2^ > 3v / 2tq' /4 / \/tt. We see that it takes longer 
for the long wavelength modes to dominate G2 than the 



order parameter G\. Because of the exponential growth, 
Pi and P2 are 0(1). We stress that these are lower 
bounds. 

If the linear equation (^7j) were valid for later times 
At > At2 then the integrals are dominated by the saddle- 
point at t' — At, to give a separation of line zeroes 
£ zero (At) of the form 



eLo(At) 



Gi(At) _ 4At 
G 2 {At) ~ ~3T 



i6 



1/2 



AAt 
~3T 



te q (t C ) 2 

(28) 



approximately independent of the cutoff I for I = O(£o)- 
Only then, because of the transfer of power to long wave- 
lengths, do line zeroes become widely separated, and 
£, zero (At) does begin to measure vortices. Thus, if the 
order parameter is large enough that it takes a long time 
(in units of t) for the field to populate its ground states 
( p5| ) we would recover the result £ zer o = 0(t; eq (tc)) as an 
order of magnitude result from (Eq) even though, a pri- 
ori, they are unrelated. This explains how experiments 
can be in agreement with the simple causal predictions. 

Whether we have time enough depends on the self- 
coupling (3, which determines when the linear approxima- 
tion fails. At the absolute latest, the correlation function 
must stop its exponential growth at At = At sp , when 
(|(?!>| 2 }, proportional to G\, satisfies (p5|). Let us suppose 
that the effect of the backreaction that stops the growth 
initially freezes in any defects. This then is our prospec- 
tive starting point for identifying and counting vortices. 

It happens that (|25|) is satisfied when 



G x {At sp ) = 7T 2 a 2 Q ^//3k B T c = tt 2 /^/i-T g /T c , (29) 

where Tq is the Ginzburg temperature introduced ear- 
lier. 

In order that G\ is dominated by the exponentially 
growing modes, so as to recover the modification of pre- 
diction (0) in the form (|2^) for identifiable vortices, the 



(30) 



condition At sp > At\, A<2 becomes, from (f29| 
(T Q /T Q )(l-T G /T c )<Cn\ 



on restoring tq, where C ~ 1. 

How strongly the inequality should be satisfied in (|3£ 
is not obvious, assuming as it does that the backreac- 
tion is effectively instantaneous, and that mean field crit- 
ical indices are valid. However, there is no way that 
the inequality can be remotely satisfied for He, when 
subjected to a slow mechanical quench, as in the Lan- 
caster experiment, for which tq/tq = O(10 10 ), since the 
Ginzburg regime is so large that (1 — Tq/T c ) — 0(1). 
That is, field growth must stop long before vortices are 
well defined. Further, since the experiment leaves the 
superfluid in the Ginzburg regime, thermal fluctuations 
will inhibit the creation of stable defects, although it may 
be that incoherent fluctuations, even if not vortices, will 
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give a signal. Adopting renormalisation improved critical 
indices cannot repair such a deficit. Whatever, there is 
no reason to expect a vortex density (^2|) . 

The situation for 3 He — B is potentially very different, 
For rapid quenches a TDLG approach is valid G^] . Sup- 
pose that (|3C]), which looks like an inequality between the 
quench rate and the equilibriation rate that permits un- 
stable modes time to dominate over short-range fluctua- 
tions, remains true, albeit with new coefficients. Firstly, 
for 3 He — B, the Ginzburg regime is very small, with 
1 - T G /T C = O{10- 8 ). The quench ends outside it. Sec- 
ondly, in generating the phase transition by nuclear reac- 
tions, rather than by mechanical expansion, the quench 
rate is increased dramatically, with tq — 10 2 — 10 4 . The 
inequality (|3p| ) is satisfied by a huge margin and we can 
understand the success of the Helsinki and Grenoble ex- 
periments. 



B. QFT: mode growth v fluctuations 

It is because the formation of defects is an early-time 
occurrence that it is, in large part, amenable to analytic 
solution. Again we revert to the mode decomposition of 
Eq.|l9|. The field becomes ordered, as before, because of 
the exponential growth of long-wavelength modes, which 
stop growing once the field has sampled the groundstates. 
What matters is the relative weight of these modes (the 
'Bragg' peak) to the fluctuating short wavelength modes 
in the decomposition Eq.|2^ at this time, since the con- 
tribution of these latter is very sensitive to the cutoff I. 
Only if their contribution to Eq.|ll| is small when field 
growth stops can a network of vortices be well-defined at 
early times, let alone have the predicted density. Since 
the peak is non-perturbatively large this requires small 
coupling, which we assume. 

However, to determine the counterpart to Eq.[30] for 
QFT is not so simple, and it helps to begin with an 
idealised problem, in which the quench is instantaneous 
(t Q = 0), and m 2 (t) flips from +e M 2 to -M 2 at t = 0. 
Eq.[n] is clearly inappropriate, a useful warning that 
the simple picture above has its limitations. However, 
once we have understood the instantaneous case the slow 
quench is not much more difficult. 

The mode equations Eq.|l^ are instantly solvable, and 
G(r, t) of Eq|2(] is simple calculable. 

Rather than introduce a cutoff at scale I, l^ 1 = A = 
0(M), as in Eq.pT|, we adopt the more brutal approach of 
cutting off modes at |k| < A in EqJ|^. Imposing the KMS 
boundary conditions at t < for an initial Boltzmann 
distribution at temperature Tq compatible with ther- 
mal mass y/eoM determines the resulting coarse-grained 
G*A(r;i) uniquely. For t < all modes are oscillatory, 
but for t > long wavelength modes with k 2 > A 2 are 
unstable, and grow exponentially. For A = 0(M), but 
A > M, say, GA(r;t) can be decomposed in an obvious 
way as G\(r;t) = G™(r) for t < 0, and 



G A (r; t) = GT(r) + G| k | <M (r; t) + G A>M>M (r; t) (31) 

for t > 0, where the second and third terms are the 
long and short wavelength fluctuations that grow from 
t = onwards. 

The equilibrium background term G™ (r) has the form 

G% l (r) = [ fke ik *C(k), (32) 

J|k|<A 

where C{k) encodes the initial conditions as before. The 
correlation length of the field is £o = M -1 and, for cutoff 
A = 0(M), the cold vortex thickness, we have line zero 
density n zero = 0(A 2 ), as in the condensed matter case. 
There is a sea of line zeroes, largely tiny loops, separated 
only by the Compton wavelength. The oscillating field 
fluctuations make them unsuitable as serious candidates 
for string since the density depends critically on the scale 
A -1 at which we view the field. 

For early positive times n zero is equally sensitive to 
A but, once Mt 3> 1, the relevant term is G|k|<Af C 7 "! *)) 
whose integral at time t is dominated by the nonpertur- 
batively large peak [Q in the power of the fluctuations 
at k around ko, where tk 2 . ~ M. Once ko <§C M we have 
the required insensitivity of the line density to the scale 
at which we coarsegrain. We find, approximately, that 

G\ kl<M (r; t) (x ^^e 2Mt exp{-r 2 M/40. (33) 

The prefactor comes from approximating C(k ) by 
To/eoM 2 . For t > the unstable modes grow until 
G| k |< M (0;^p) = 0{4>l) = 0{M 2 /X), determining the 
spinodal time t sp as Mt sp = 0(ln(l/A)). The long- 
distance behaviour of G™(r) = 0(e~ Mr ), with its cor- 
relation length £ = M -1 is a shortlived relic of the initial 
thermal conditions. After a time t r — 0(M^ 1 ) < t sp , 
it is rapidly supplanted by the behaviour of the expand- 
ing long wavelength modes G| k | <M (r;£)/G| k | <M (0;t) w 
exp{-r 2 /£ 2 (i)} where £, 2 (t) = tt/M « 4/fc 2 . With 
Gi k i < ^f(0;t S j,) = 0(A _1 ) non-perturbatively large and 
GX"(0) and G A> | k | >M (0;i sp ) of order A -1 / 2 , 

n Z ero(t sp ) « ^ a J (1/A) [l + 0(X 1/2 ln(l/A))]. (34) 

The relative error term in the brackets is, in large part, 
a measure of the fluctuation vortices that we mentioned 
earlier and is a measure of the stability of this density 
to changes in the coarse-graining scale A = 0(M). For 
weak coupling 0(A 1/2 ln(l/A)) < 1. 

There are other, less direct, ways to understand why 
the strings only become well-defined once Ga/(0;£) is 
nonperturbatively large, even though the density is in- 
dependent of its magnitude. Instead of a field basis, we 
can work in a particle basis and measure the particle 
production as the transition proceeds. Whether we ex- 
pand with respect to the original Fock vacuum or with 
respect to the adiabatic vacuum state, the presence of a 



8 



non-perturbatively large peak in k 2 G(k\t) signals a non- 
perturbatively large occupation number N^ oc 1/A of 
particles at the same wavenumber fco E3]. With n zero 
of ([54]) of order k 2 this shows that the long wavelength 
modes can now begin to be treated classically. 

From a slightly different viewpoint, the Wigner func- 
tional only peaks about the classical phase-space trajec- 
tory once the power is non-perturbatively large [ |l9|j20| 
from time t sp onwards. More crudely, the diagonal den- 
sity matrix elements (field probabilities) are only then 
significantly non-zero for non-perturbatively large field 
configurations <f> oc A -1 / 2 like vortices. 

So far we have done no more than estimate the den- 
sity of coarse-grained line zeroes for a free roll at the 
time t sp at which nothing has frozen in. For this global 
0(2) theory the damping of domain growth occurs by the 
self-interaction effectively forcing the negative m 2 (t) to 
vanish so as to produce Goldstone particles. This initial 
slowdown leads to no qualitative change in the vortex 
density. 

All the results above were for the instantaneous 
quench. We now return to the original problem of slower 
quenches, with e(t) as in Eq.|lJ, in which the symmetry- 
breaking begins at relative time At = 0.. For a free 
roll, the exponentially growing modes that appear when 
At > tZ = tQk 2 /M 2 lead to the approximate solution 



G(r;At) 



T 



( M 



M\m(At)\ Vv/AH^ 



■V 2 4MAt 3 / 2 



-r 2 /e(At) 



(35) 



where £ 2 (Ai) = 2-J AUq /M . The provisional relative 
freeze-in time At sp is then, for MtQ < (1/A), 

MAt sp ~ (A« Q ) 1 / 3 (ln(l/A)) 2 / 3 ~ Aft c (ln(l/A)) 2 / 3 . 

(36) 

This is greater than Mtc, but not by a large multiple 
unless we have a superweak theory. More detail is given 
in@. 

At this qualitative level the correlation length at the 
spinodal time is 



M 2 e(t sp ) c (Mt Q ) 2 / 3 (ln(l/A)) 1 / 3 . 



(37) 



The effect of the other modes is larger than for the in- 
stantaneous quench, giving, at t = t sp 



M 2 



tt(Mtq) 2 / 3 



(ln(l/A))- 1 / 3 [l + S] 



(38) 



The error term E = 0(A 1 / 2 (Mi Q ) 4 / 3 (ln(l/A))- 1 / 3 ) is 
due to fluctuation vortices. In mimicry of EqJT| it is 
helpful to rewrite Eq.Bq as 



2/3- 



(ln(l/A))- 1 /3[i + ^]. (39) 



in terms of the scales tq = £0 = M~ x . The first term in 
Eq.|39| is the Kibble estimate of Eq.|lT|, the second is the 
small multiplying factor, rather like that in Eq.|2^, that 
yet again shows that estimate can be correct, but for com- 
pletely different reasons. The third term shows when it 
can be correct, since E is also a measure of the sensitivity 
of n zero to the scale at which it is measured. We note that 
if we take Mtq = ln(l/A) we recover the instantaneous 
results qualitatively. Only for larger tq will the quench 
give different results. The condition E 2 -C 1, necessary 
for a vortex network to be defined, is then guaranteed if 



(t q /t q ) 2 (1-T g /T c )<1, 



(40) 



on using the relation (1 — T G /T c ) — 0(A). This is the 
QFT counterpart to Eq.||. 

The density of Eqj39| is smaller than that of Eq.|3^ by 
the factor (ln(l/A)) 1 ^which is, in principle, a large num- 
ber even if in practice it is only a small multiple. The 
easiest way to enforce E <C 1 and Mtq > ln(l/A) is to 
take Mtq = ln(l/A) Q , for a > 1. The effect in Eq.|| is 
merely to renormalise the critical index. Of course, the 
Kibble prediction Eq.[ll] was only an estimate. Although 
it is good qualitatively it is misleading when consider- 
ing definition of the network since a simple calculations 
shows that, at time tc, the string density is still XSto- 
tally sensitive to the definition of coarse-graining. 

Finally, suppose that this approach is relevant to the 
local strings of a strong Type-II U(l) theory for the 
early universe, in which the time-temperature relation- 
ship tT 2 = YMpi is valid, where we take Y = C^ICT 1 ) 
in the GUT era. If G is Newton's constant and \i 
the classical string tension then, following ||, MtQ ~ 
10 _1 A 1//2 (G/Lt) -1 / 2 . The dimensionless quantity G[i ~ 
1CP 6 — 10~ 7 is the small parameter of cosmic string the- 
ory. A value A — 1CT 2 gives Mt Q ~ (Mt sp ) a ,a ~ 2, 
once factors of n, etc. are taken into account, rather than 
MtQ ~ 1/A, and the density of EqJH may be relevant. 

We have not attempted to show what happens directly 
after the formation of defects for those cases when they 
are well-defined by the time the initial population of the 
ground states has occurred. We have only attempted to 
test the Kibble/Zurek scenarios for defect formation. For 
strongly-damped systems like condensed matter systems 
the linear approximation seems to be valid for longer 
times than might have been guessed. That would explain 
why the final population is simply related to the initial 
population. For QFT it is much less obvious that this is 
the case. That is, even when Kibble is correct, the later 
evolution may make the estimate irrelevant. To take one 
extreme example. If we suppose that the field evolution 
can be approximated by self- consistent linearisation then, 
very rapidly, the scattering between field modes, absent 
in the linear approximation here, will be enough to re- 
distribute power back to short wavelengths. The density 
of line zeroes then begins to increase rapidly. If they do 
represent strings their density is likely to be independent 
of tqH. 
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VI. CONCLUSIONS 

The Kibble/Zurek scenario for the onset of a continu- 
ous transition, either in QFT or condensed matter, sug- 
gests that causality determines the correlation length of 
the fields at the time that they freeze in. Our simple 
calculations show that this is qualitatively correct. 

A further step in their arguments was to identify this 
calculable correlation length with the length of defect 
separation (cosmic strings or vortices) at the time at 
which such defects were produced. Whereas this is in 
agreement with vortex experiments on 3 He it is not the 
case for 4 iJe. In order to understand this we have iden- 
tified defect density with field-zero density, once this is 
effectively independent of the scale at which the zeroes 
are counted. 

The separation of zeroes is determined from the short- 
distance field correlations and the correlation length (by 
definition) from the long-distance field correlation func- 
tion. Nonetheless, the field correlation length at the time 
it freezes in and the zero separation length when the ini- 
tial relaxation of the fields to their ground states is com- 
plete are qualitatively the same, provided thermal fluctu- 
ations are not too strong for a well-defined vortex net- 
work to be established. Qualitatively, this requires (in 
the terminology of the text) 



(T Q /T y(i-T G /T c )<i, 



(41) 



where (in mean field) 7 = 1 for condensed matter, and 
7 = 2 for QFT. 

This linking of quench rate to the Ginzburg regime 
explains why 3 He and 4 He experiments are so different, 
and why we expect the simple predictions for vortex den- 
sities not to be verified in the latter. A constraint like 
Eq.[ll] is not strong enough to preclude initial defect for- 
mation in the early universe. 



ACKNOWLEDGEMENTS 



[is; 

[19 

[20 

[21 
[22 



[4] W.H. Zurek, Nature 382, 297 (1996), 

P. Laguna and W.H. Zurek, Phys Rev. Lett. 78 2519 
(1997). A. Yates and W.H. Zurek, |hep-ph/980i223 . 

[5] M.M. Salomaa and G.E. Volovik, Rev. Mod. Phys. 59, 
533 (1987) 

[6] V.M.H. Ruutu et al., Nature 382, 334 (1996). 
[7] C. Bauerle et al, Nature 382, 332 (1996). 
[8] P.C. Hendry et al, Nature 368, 315 (1994). 
[9] P.V.E McClintoch et al, to be published in Phys. Rev. 
Letters. 

[10] G. Karra and R. J. Rivers, to be published in Phys. Rev. 

Letters. 

[11] M. Bowick and A. Momen, preprint }iep-ph/9803284| , to 

be published. 

[12] V.N. Popov Functional integrals and collective excita- 
tions, Cambridge University Press (1997), Chapter 8. 

[13] T. Vachasparti and A. Vilenkin, Phys. Rev. D30, 2036 
(1984), A. Vilenkin and R.J. Scherrer, Phy s. Rev. D5 647 
(1997), see also preprint friep-ph/970949 > fj 

[14] M. Gleiser and H-R. Muller, preprint tiep-lat/970400^ , 
to be published in Physics Letters B. 

[15] B.I. Halperin, published in Physics of Defects, proceed- 
ings of Les Houches, Session XXXV 1980 NATO ASI, ed- 
itors Balian, Kleman and Poirier (North-Holland Press, 
1981) p.816. 

[16] F. Liu and G.F. Mazenko, Phys. Rev. B46, 5963 ( 1992). 
[17] Yu. M. Bunko v and O.D. Timofeevskaya, preprint cond- 
mat/9706004 , Phys. Rev. Letters, to be published. 
D. Boyanovsky, Da-Shin Lee and A. Singh, Phys. Rev. 



D48, 800 (1993), 

D. Boyanovsky, H.J. de Vega and R. Holman, Phys. Rev. 
D49, 2769 (1994). 

A. Guth and S-Y. Pi, Phys. Rev. D32, 1899 (1985). O. 
Eboli, R. Jackiw and S-Y. Pi, Phys. Rev. D37, 3557 
(1988). 

S. Mrowczynski and B. Muller, Phys. Rev. D50, 7542 
(1994). 

G. Karra and R.J. Rivers, Physics Letters B414, 28 
(1997) 

H. J. de Vega et al., private communication. 



I particularly thank Glykeria Karra, with whom most 
of this work was done. I also would like to thank Alasdair 
Gill, Tom Kibble, Wojciech Zurek and Yuriy Bunkov for 
fruitful discussions. This work is the result of a network 
supported by the European Science Foundation. 



[1] T.W.B. Kibble, J. Phys. A9, 1387 (1976). 

[2] T.W.B. Kibble, in Common Trends in Particle and Con- 
densed Matter Physics, Physics Reports 67, 183 (1980). 

[3] W.H. Zurek, Nature 317, 505 (1985), Acta Physica 
PolomcaB24, 1301 (1993). See also W.H. Zurek, Physics 
Reports 276, Number 4, November 1996 



10 



